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Abstract 

Using multi-dimensional PIC simulations we show that the magnetic undulator-type field of the 
plasma magnetostatic mode is indeed produced by the interaction of a laser pulse with a relativistic 
ionization front, as predicted by linear theory for a cold plasma. When the front with this mag- 
netostatic mode is followed by a relativistic electron beam, the interaction of the beam with this 
magnetic field, produces FEL-type synchrotron radiation, providing a direct signature of the mag- 
netostatic mode. The possibility of generating readily detectable ultrashort wavelength radiation 
using this mode, by employing state-of-the-art laser systems, is demonstrated, thus opening the 
way towards experimental observation of the hitherto unseen magnetostatic mode and the use of 
this plasma FEL mechanism to provide a source of high-brilliance ultrashort wavelength radiation. 
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The advent of high power lasers has allowed for the exploration of the extraordinary 
ability of plasmas to sustain strong electric pQ and magnetic [2j fields. Electrostatic plasma 
waves [3], driven by intense laser pulses have well-defined electric structures that can be 
used to accelerate electron beams up to the GeV energy jIHTj. Laser-plasma interactions 
are also known to provide ultra-strong (> 100 T) magnetic fields [2j [8]. However, effective 
methods for generating controlled periodic magnetostatic structures in plasmas are yet to 
be demonstrated. 

In a recent study, it was shown that an overdense ionization front moving nearly at c 
can frequency upshift microwave radiation confined in a cavity [9]. This study, and earlier 
theoretical work [10], has led to the understanding of how electromagnetic boundary condi- 
tions should be handled at a plane dielectric discontinuity moving normal to the plane. A 
zero-frequency, zero-electric-field, purely Ampere-law magnetostatic mode (here called the 
MS mode) was recognized as being important. This mode has hitherto not been seen and 
is ordinarily very difficult to excite. With a plasma front boundary which is moving at a 
normal velocity Vf ) the general phase continuity condition is that the quantity (uj — k • Vf) 
does not change when changing from one medium to the other in uj — k space [10]. With 
a fast-moving front and an incident EM (ElectroMagnetic) wave whose frequency uji is well 
below the plasma frequency (uj p ) behind the front, the stationary MS mode becomes easy 
to set up, and the MS mode (which has the same magnetic polarization as the incident EM 
wave) becomes the dominant repository for the incident energy [lOj. (The value of uj p is 
usually that obtained after rapid and complete ionization by a high frequency (uj^ 3> uJ p ) 
ionization front driver laser pulse of the gas behind the ionizing pulse.) The MS mode will 
have a normal wavevector kMS of (1 + v f/ c ) times the incident wavevector (in the case of 
interest this factor is very nearly 2). Since the MS mode does not propagate and cannot 
exist outside the plasma, it has not yet been detected directly. 

For a plane polarized incident EM wave the MS mode in the plasma has the magnetic 
field geometry of standard magnetic undulators, such as those used for FELs (Free Electron 
Lasers), and from this arises the key concept of this letter in addressing the MS mode 
detection problem. One begins by creating a usefully strong MS undulator-type B-field 
using the sub-plasma-frequency laser wave incident on the counter-propagating ionization 
front created by a short ionization-front laser drive pulse incident from the opposite side 
of the plasma. The key new component here is to supply also a highly relativistic electron 
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beam pulse to pass through the stationary MS undulator magnetic field almost immediately 
after it is produced. This is done in order to produce a pulse of synchrotron radiation more 
or less in the electron beam direction, thereby confirming the existence of the MS mode. 

This MS synchrotron radiation will provide both a direct diagnostic of the MS spatial 
period (from the frequency spectrum peak) and a reasonable estimate of its strength (from 
the magnitude of the synchrotron signal). Going beyond this basic verification of the MS 
mode physics brings in the complementary notion of using the concept to produce a powerful 
undulator magnetic field of much shorter spatial period than is possible to obtain with 
conventional magnets. Applications might range from condensed matter physics [HJ CE2] 
to novel radiations sources [13]. Showing how this would work in practice is the purpose 
of the rest of this paper, as exemplified by the simulation results and the further detailed 
discussions. 

The efficient generation of the MS mode requires a sharp ionization front, i.e. L g ki <C 1, 
where L g is the gradient scale of the ionization front and ki is the wavenumber of the incident 
light [10]. Steep ionization fronts can be generated by the propagation of a short intense 
laser pulse in a gas jet via tunneling ionization. The relativistic factor of a laser-driven 
ionization front, in the linear regime, is given by 7/ = Ud/uJp- For plasmas transparent to 
the ionizing light, uoa 3> and for uj^ 3> ujp/^jj the amplitude of the MS mode [10J can 
be written as 



is only limited by the threshold for ionization of the background gas. For Hydrogen, for 
instance, the threshold intensity for tunnel ionization [14] is ~ 10 13 W/cm 2 for mid-infrared 
light, allowing for the generation of magnetic-field amplitudes on the order of 100 T. 

Present state-of-the-art CO2 laser systems [15J can deliver pulses with energy q = 170 
J and = 10.6 /xm, thus allowing for the generation of a MS amplitude of 100 T, a 
wavelength of 5.3 /xm, and a controllable width and length (e.g. with spot size of 400 /xm 
and length ~ 10 cm), with the MS mode- width determined by the smaller of the spot size 
of the incident pulse, Wi, and the width of the ionization front (approximately the spot size 
of the driver pulse, W^), and the MS length is typically determined by the incident pulse 
duration, 7$. Figs. [T] and [2] illustrate the features of the MS mode generated by the collision 
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where \ = 2irc/u>i is the wavelength of the incident radiation and 7j its intensity, which 
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of a CO2 laser pulse with a relativistic ionization front in two-dimensional (2D) and three- 
dimensional (3D) OSIRIS 2.0 [IS] PIC simulations. The ionization front is self-consistently 
generated by tunnel ionization due to a short, linearly polarized laser pulse. The CO2 laser 
pulse is launched in counter-propagation with the ionization front. We have also launched 
an electron beam trailing the driver pulse in order to probe the generated MS mode. We 
can observe that behind the ionization front a modulated magnetic field is induced, with a 
wavelength of 5.3 fim and amplitudes of 54 T, for n = 5 x 10 19 cm -3 (Fig. [I]), and 93 T, 
for n = 5 x 10 20 cm -3 (Fig. [2]), in very good agreement with Eq. ([I]). Moreover, there is 
almost no transmission of the CO2 pulse in that region. The MS mode (with wavevector 
2kj) is stationary in space behind the ionization front, which itself moves to the right at ~ c, 
while the incident field (with wavevector kj) moves to the left at c (Fig. [TJd) . The transverse 
momentum of the beam electrons, p x: shown in Fig. [2| illustrating the wiggling motion in 
the MS mode, is in very good agreement with the theoretical prediction |p x | max = K [TT] . 
where K = 0.934£> u [T]A u [cm] is the dimensionless undulator strength parameter, B u is the 
amplitude of the undulator magnetic field, and \ u is the wavelength of the undulator. The 
presence of the trailing e-beam does not affect the MS structure, indicating that the MS 
mode can be used as a magnetic undulator. 

As the beam electrons wiggle in the undulator structure they emit radiation with wave- 
length A r = \ U {1 + K 2 /2)/2^l [T7j, where 75 is the relativistic factor of the beam electrons. 
The short wavelength (X u ~ 5 /im) and high amplitude (B u <~ 100 T) laser driven MS struc- 
ture allows for the generation of radiation with features that can provide an evidence of the 
MS mode (energy determined by the properties of the e-beam and MS structure, narrow 
spectrum, and high brilliance). The synchrotron radiation emitted by the e-beam traversing 
the MS mode, and therefore the MS mode itself, can be unambiguously detected through a 
straightforward null test whereby the e-beam/incident laser pulse are present or not in the 
experimental set-up. 

The total power radiated by an e-beam traversing the MS undulator-like field is given 
by Q^P/e, where P = (l/3)e 4 £> 2 7 2 /m 2 c 3 [17J is the power radiated by a single electron 
in the undulator structure, yielding P r [kW]~0.71 (B u [T}/60) 2 (e 6 [MeV]/100) 2 (Q 6 [nC]/0.1). 
The average radiated energy is P r r p /2 ) where r p is the propagation time of the e-beam in 
the MS structure. Since the residence time of the e-beam in the MS structure is half the 
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duration of the incident wave, 7$, the total radiated energy can be written as 

er[ „j] * 28^49 (*m 2 f«V 

L J V 60 y v ioo / o.i 160 

The energy of the radiated photons is 
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e7 [keV]. 18.72^ («Y, (3) 
7L J Xilfim] V 100 J K J 

which depends only on the energy of the e-beam and on the MS wavelength. We note that 
while the MS magnetic field can be strong, K is still small due to the short wavelenght 
of the MS undulator, and, therefore, the photon energy is well approximated by 47 2 times 
the incident laser photon energy. The total number of photons, N 1 = e r /e 7 , is given by 
7V 7 [10 7 ] ~ 0.95(5 u [T]/60) 2 (A,[/xm]/10.6)(g 6 [nC]/0.1)(r,[ps]/160). For typical parameters, 
that can be reached in an experiment, 10 7 photons are emitted in a half cone angle of 1/75, 
and can be detected far away from the plasma, using an x-ray CCD camera in a single 
photon mode (where as few as 3-5 photons can be detected in the energy range of interest) 
and suitable spectral filtering. It is also important to point out that the e-beam propagates 
along a large number of periods, N U) in the undulator-like structure. Thus, for a low energy 
spread e-beam, Ae^/e^ <C 1, the energy spread of the emitted radiation, Ae 7 /e 7 , is narrow, 
since it is determined by the largest of 1/N U or Ae^/e^. The dimensionless peak brilliance 
B 1 of the radiated photon beam, normalized to 1 photon/s/mrad 2 /mm 2 /0.1% bandwidth, 
is5 7 [10 28 ] ~ 1.52( J B u [T]/60) 2 (r i [ps]/160) 2 (3/r 6 [^m]) 2 (0.1/^[mrad]) 2 (g 6 [nC]/0.1)(10/r 6 [fs]), 
where is the radius of the e-beam and 9 its angular divergence. The propagation of an 
e-beam with 65 = 1 GeV through the plasma undulator induced by a CO2 laser pulse, would 
provide gamma-ray pulses with e 7 ~ 2 MeV and a peak brilliance ranging from B 1 ~ 5 x 10 25 , 
for typical laser wakefield accelerator (LWFA) e-beams (Q& = 30 pC, 9 = 1 mrad, r& = 10 fs, 
r\) = 3 /im), to B 1 ~ 10 30 for e-beams from typical conventional linear accelerators (Q& = 2 
nC, 9 = 1 /xrad, 75 100 fs, 60 /im). 

A proof-of-principle experiment, demonstrating and measuring the MS mode can be 
performed using the widely available 1 J class Ti:Sapphire laser systems, whose laser pulses 
can efficiently generate both the relativistic ionization front and LWFA e-beams with 65 = 50 
MeV and = 30 pC [5]. Using an improved two-stage Raman shifting scheme with > 0.2% 
efficiency [18], a fraction of the laser energy (0.2 J) can be used to produce pulses with 
Ti = 4 ps and = 5 /im, which then collide with the ionization front at incident intensities 



Ii ~ 3 x 10 13 W/cm 2 , Wi = 10 /xm, inducing the MS mode. With these parameters, 
N 7 ~ 10 5 , peaked at e 7 ~ 10 keV, which would allow for a clear evidence of the MS mode. 

Since the e-beam probing the MS mode propagates in a plasma, it is important to assess 
the possible beam quality degradation, in terms of energy spread and emittance, due to 
beam-plasma and laser-plasma instabilities. The relevant beam-plasma instabilities are the 
two-stream instability (longitudinal) and the Weibel instability (transverse). The two-stream 
instability can be avoided using short e-beams, L&[mm] < (e&[MeV]/100) (n5[10 17 cm -3 ]) -1 / 3 
(n e [5 x lO^cm" 3 ])" 1 / 6 [19], where rib is the electron density of the beam. The Weibel insta- 
bility is strongly suppressed for e-beams narrower than the plasma wavelength X p = 2ttc/uj p 
[20] . or for wider low-density e-beams with n b /n e < 2 x 10- 3 (e 6 [MeV]/100)(%rad]/10) 2 [21]. 
In what concerns laser-plasma instabilities, laser filament ation is avoided for laser powers 
smaller than the critical power for self- focusing [22] . and stimulated Raman scattering [23] 
will not be driven for ultrashort laser pulses, which are required in order to generate a 
steep ionization front. Laser wakefields [H [24], which can also have deleterious effects on 
the e-beam quality, are negligible for a laser normalized vector potential a ^ <C 1, since 
5n e /n e0 <C 1. 

Ionization front irregularities can also contribute to a degradation of the MS mode and 
consequently of the emitted radiation. Erosion of the head of the laser pulse can broaden 
the initially sharp ionization front and, therefore, the laser should contain enough energy 
to generate a stable ionization front during the interaction process. The energy depletion 
length, Ldpi, due to ionization can be roughly estimated by equating the laser energy to 
the ionization energy of a single atom, Wi on: times the number of atoms ionized, n^a^T^ 
(where aa is the laser cross section), yielding Ld p i[cm]~ 10/^(2 x 10 16 Wcm _2 ]r^[50fs]/(n [5 x 
10 19 cm _3 ]W^ on [13.6eV]). For typical experimental parameters, Ld p i exceeds several centime- 
ters. Phase irregularities across the ionization front can also lead to a local variation of 
the MS wavelength, increasing the energy spread of the emitted radiation. However, near 
the axis this energy spread induced by phase irregularities should be very small (the energy 
spread is proportional to 1— cos[0], where is the variation in angle of the tilted wavefront), 
and therefore, for e-beams that are narrower than the ionization front (driver pulse) this 
should not significantly affect the brilliance of the emitted radiation. 2D simulations of the 
ionization front propagation in lengths in excess of 1 cm have confirmed this, and will be 
presented elsewhere. 
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We have performed ID PIC simulations fully resolving the radiation generated due to 
the interaction of the relativistic e-beam with the induced MS mode, with a numerical 
subtraction technique to isolate the synchrotron radiation of the e-beam; by running two 
identical simulations, with and without the e-beam, the radiated fields are obtained from the 
subtraction of the results of the two simulations. Fig. [3] illustrates the temporal evolution 
of the energy radiated by the e-beam, showing the generation of radiation with e 1 = 45.2 
eV. The induced MS amplitude varies from 93 T (when the e-beam starts to interact with 
the structure) to 69 T (at the end of the interaction), due to the slow increase of the rise 
time of the ionization front as the driver pulse propagates [25]. The total radiated energy 
and the radiation spectrum are in good agreement with Eqs. ^ and 

In conclusion, we have proposed a configuration that allows for the detection and char- 
acterization of the hitherto unseen plasma MS mode. It is induced by the collision of light 
pulses with relativistic ionization fronts, and detected by using it as a plasma undulator. 
PIC simulations show the generation of periodic, short wavelength magnetic fields associated 
with the excitation of the MS mode, and illustrate the generation of ultrashort wavelength 
synchrotron radiation by the propagation of an electron beam through this magnetic struc- 
ture. Our results open the way towards the first experimental evidence of the MS mode 
from the measurements and analysis of the emitted radiation. Once this is accomplished, 
the possibility of using the plasma MS mode with an electron beam to produce an effective 
and useful source of ultrashort wavelength radiation can then be pursued. 
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FIG. 1: (color online) a) Snapshot of the MS mode generated by a CO2 laser pulse (Aj = 10.6 /im, 
ei = 8 mJ, Wi = 165 /im, = 1.25 ps, /j~3x 10 13 W/cm 2 ) in 2D PIC simulations. A short pulse 
(Ad 0.4 /im, 6^ = 2 mJ, 50 /im, 1^ — 2 x 10 15 W/cm 2 , 7^ 50 fs) is launched from the left 

hand side of the simulation box through an Hydrogen gas jet (no = 5 x 10 19 cm -3 ) and is trailed 
by an e-beam (energy 65 = 100 MeV, charge — 10 pC, duration 75 = 25 fs, radius — 5 /im, 
and density = 10 17 cm -3 ). The CO2 laser pulse is launched from the right hand side of the 
simulation box. The simulation box size is 450 x 450 /im 2 , with 2 15 x 2 11 cells and 9 particles 
per cell (roughly 600 million particles were used), b) Time streak of a longitudinal lineout of the 
magnetic field at the center. The arrows and the dashed line indicate the position of the front. 




FIG. 2: (color online) 3D snapshot of the MS mode. Beam electrons are colored by transverse 
momentum and the B-field of the CO2 laser is truncated at the ionization front position. The laser 
and e-beam parameters are identical to those in Fig. [I] for a scaled simulation box with a width of 
180 /im and a length of 450 /im. 
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FIG. 3: (color online) Energy radiated by the propagation of an e-beam (e& = 5 MeV, 75 = 10.5 fs, 
715 = 8.2 x 10 16 cm -3 ) through the MS mode (solid), generated by the collision of a CO2 laser pulse 
with a relativistic ionization front, in ID PIC simulations. The laser parameters (ionization- front 
driver and incident CO2 laser) are identical to those in Fig. [l]and the Hydrogen gas jet has a neutral 
density no = 5 x 10 20 cm -3 . The simulation box is 1.5 mm long, with 2 19 cells and 100 particles 
per cell (total ~ 52 million particles). The theoretical radiated energy, Eq. ([2]), is plotted for the 
initial (dashed) and final (dot-dashed) B-field amplitude. After 1.8 ps of propagation (N u = 100), 
the radiation spectrum exhibits a narrow energy spread Ae 7 /e 7 = 0.012 (shown in the inset). 
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